Irreversible shape distortions by flux-pinning-induced magnetostriction in hard 

type-II superconductors 
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Exact analytical results are obtained for the flux-pinning-induced magnetostriction in long cylinders 
of type-II superconductors placed in a parallel magnetic field. New effects are found regarding 
the irreversible deformation of the cross-sectional area: Whereas the magnetostriction of a circular 
cylinder is shape conserving, it is shown that a square cross-section deforms with considerable 
distortion. During a field cycle both concave, convex, and even more complicated distortions are 
predicted. The strong implications for the interpretation of dilatometric observations are pointed 
out. The main results obtained for both the circular and square geometry are valid for any critical- 
state model, jc = jc{B), and we use the Bean model for illustration and discussion. 

PACS numbers: 74.25. Ha, 74.60. Ge, 74.60. Jg 



It was discovered by Ikuda and coworkers^ that high 
temperature superconductors (HTSCs), when placed in 
a magnetic field, can respond by a giant magnetostric- 
tion. Relative changes in sample size as large as 10~'* 
were observed in single crystals of Bi2Sr2CaCu208. To 
explain this deformation, which is about two orders of 
magnitude largely than values reported for conventional 
superconductors p they considered the stress exerted on 
the crystal by pinned vortices. Treating the crystal as 
an infinite slab of half- width, w, placed in a parallel ap- 
plied field, Ba/ fJ-o, the pinning-induced magnetostriction 
in the transverse direction was derived from the critical- 
state model to be 



Aw/w = (2cfiowy 



[B\x)~Bl]dx, 



(1) 



where c is the elastic constant and B{x) the local in- 
duction. As evident from eq.(l|), and also pointed out 
immediately by Ikuta et al., the pinning-induced mag- 
netostriction is intimately related to the magnetization, 
and consequently should display similar irreversible char- 
acteristics. Following the discovery, the giant magne- 
tostriction was observed also in single crystals of other 
HTSC materialsja thus establishing the effect as being a 
general phenomenon which, like the magnetization, can 
provide important information about the pinning prop- 
erties. 

Magnetostriction, as dilatations in general, are most 
frequently measured by displacive sensors like e.g. a ca- 
pacitance dilatometer cell. Such sensors provide an elec- 
trical signal which normally represents_the variation in 
the end-to-end dimensions of a samplcQ Although very 
high sensitivities (10^^^ m) can be achieved, the cells 
are usually designed so that distortions in shape are not 
accounted for, or they may even influence the measure- 
ments in an uncontrolled way. Up to now, shape dis- 
tortions in connection with the pinning-induced magne- 
tostriction were never considered, and all data analysis 



has been based on the relation eq. (0) valid for plane strain 
deformations only. 

In this letter we show that the pinning-induced stress 
indeed leads to significant shape distortions even in the 
simple geometry of an infinitely long isotropic cylin- 
der with square cross-section. For comparison we also 
present the exact solution of the magneto-elastic prob- 
lem for the circular cylinder where shape is conserved. 

From electrodynamics the force which drives the vor- 
ticies into the material is given by j x B per unit volume. 
For HTSCs in the mixed state under practical conditions 
the current density, j, is usually set equal to (V x B)//io 
as the|-thermodynamic field can be well approximated by 
'B/fiQU In the critical-state model one has |j| = jc, and 
a non-uniform distribution of pinned fiux transmits onto 
the crystal a volume force density equal to f = jc x B. 
In the two geometries under consideration the symmetry 
of the force field is quite different, see Fig.l. In contrast 
to the simple radial distribution shown in (a), the square 
cylinder in (b) is divided into four regions each containing 
distributed forces pointing the same direction. Moreover, 
in both cases the cylinder generally has additional divi- 
sions into shells where the forces alternate in sign - or 
they can even vanish - all depending on the magnetic 
prehistory. In all such states the force exerted on the 
material per unit volume is given by 



f = (V X B) X B//io = -(2/io)"^ VB^ 



(2) 



The last form, expressing a purely magnetic free energy 
of the vortex lattice, is used in the present calculations. 
In the analysis the strains are assumed well below the 
fracture limit allowing linear elasticity theory, and hence 
the principle of superposition, to be applicable. 

In the circular cylinder the rotational symmetry im- 
plies that the deformation is described by a radial dis- 
placement field, Ur{r). The non- vanishing components of 
the strain are e^ = uj,(r) and eg — Ur/r, which are related 



to the stress components, ag and Ur by Ecr = (Jr — i^o'g 
and Eeg = ae — vdr- Here E and v is the Young's mod- 
ulus and Poisson's ratio, respectively. 
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FIG. 1. Flux penetration and current flow pattern in a cir- 
cular (a) and square (b) cylinder in the critical-state model. 
The indicated magnetized states contain a flux front, illus- 
trating that the cylinders generally are divided in shells with 
different volume forces. 



The condition of static equilibrium is that (dar/dr) + 
(cr,. — iye)/r + / = 0, which in terms of the displacement 
field is expressed as 



u^ 



f{r) = . 



r^ E 

Integration twice, and using that Ur{0) — 0, gives 

l-J/2 1 



Ur-{r) = -- 



E 



r'F{r')dr' + Cr , 



(3) 



(4) 







where F{r) = J^ fir')dr' = [^^(q) _ B^{r)]/2fio. With 
the constant C determined by the free surface condition 
(Tr{r — R) = 0, the result becomes 



Ur{r) 







r'B\r')dr'\ , 



r'B^{r')dr' 



(5) 



which is the exact solution from which a complete stress- 
strain picture can be expressed in terms of the flux den- 
sity B{r). E.g., the radial stress is given by 



"'M-i^'«-«.-+^ 



r'B^{r')dr' 



^ r r'B\r')dr' \ 



(6) 



The eqs.(^ and (|sb| generalizes expressions for Ur and 
Or derived recentlyQ only for the Bean model, jc{B) — 
constant. We emphasize that our results are valid for any 
critical-state model. 

For the external dilatation, AR/R = Ur{R)/R, one 
obtains from eq.(||) 



AR 
R 



EfioR'^ 



r[B-'ir)~Bl]dr. 



(7) 



This formula for the pinning-induced magnetostriction 
of a circular cylinder, first found by Johansen et al.JII, is 
very similar to the expression in eq.(0) for the slab. The 
extra factor r in the integral is a simple modification, 
which one finds analogously in the expressions for the 
critical-state magnetization. The factor 1 — i/ reflects the 
stronger constraint on free expansion for the cylinder as 
compared to the slab. 
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FIG. 2. Magnetostriction, ARjR x (1 - v){^iojoR)'^ / E^o, 
as function of applied field for a circular cylinder with jc = 
Jo exp(— |_B|/_Bo). Curve (1); Bo = oo i.e. the Bean model 
and (2); Bo = /iojo-R/0.3. The field axis is normalized by the 
full penetration field Bp = Bo ln(l + nojoR/Bo). 

From eq.(0) analytical expressions for the various 
branches of magnetostriction hysteresis loops are read- 
ily derived. Figure 2 shows how AR/R depends on the 
applied field for two commonly used critical-state mod- 
els. As Ba runs through a complete cycle AR/R displays 
a huge hysteresis loop where in the Bean model case both 
the main ascending and descending branches are linear 
provided a sufficiently high maximum field. Among other 
features of the Bean model results is that the maximum 
remanent magnetostriction is given by 



(Ai?/i?)„m = (1 - iy){12E)-^ fioJc R" 



(8) 



Since the remanent state has no contribution from re- 
versible components, the above relation provides a means 
to infere jc from magnetostriction measurements. In the 
field range where the main branches are linear the verti- 
cal width of the hysteresis loop equals 

{AR/R)^ ~ {AR/R)^ =2(1- y){'iEy^ Rjc Ba ■ (9) 

Note in particular that if the hysteresis width is normal- 
ized by {AR/R)^cai the value of jc can be found from the 



slope in a plot versus Ba without knowledge of the elastic 
constants. The normalized width is simply ^Ba/ ^ii^ijcRi 
a result analogous to the extensively used Bean model 
relation between jc and the width of the magnetization 
loop. Conversely will measurements of magnetostriction 
allow determination of the elastic parameter (1 — v)/ E if 
jc is known from other measurements. 

The more curved hysteresis loop resulting from the 
exponential model, jc = jo exp[— |i3|/i3o], graph (2) in 
Fig. 2, is usually a better descrupiion of observed magne- 
tostrictive behavior in HTSCsE^eI. An extensive analysis 
with various jc(i?)-functions is subject of a forthcoming 
paper, where also the rather lengthy expressions for the 
different branches in Fig. 2, is presented. 

A final remark on the circular case is that the inte- 
gral in eq.(0) also can be written as /q r-^f{r)dr, i.e., 
the AR/R is given as the second moment of the vol- 
ume force distribution just as the jc-loops determine the 
magnetization. One could be tempted to generalize the 
relation to other cylindrical geometries, and thus allow 
for a simple punified treatment as done in magnetization 
calculations .a However, such generalization would be in- 
correct, as will become evident as we now turn to the 
square cross-section case. 

A different approach is chosen to analyze the magneto- 
elastic problem of the square cylinder. Within the linear 
elasticity approximation one can separate the treatment 
into an infinite number of deformations produced by in- 
finitessimal square loops of the force field (Fig.l b). The 
final result is obtained by superimposing the elastic re- 
sponse of the whole cylinder to each force loop. For this 
to be a useful approach one needs a way of adding these 
deformations in a coherent manner. We first show that 
this indeed is feasable for the square cylinder. 

Consider the effect of a square force loop of width dx 
with sides a distance x from the origin (Fig. 3). The ma- 
terial enclosed by the loop will experience only a nor- 
mal stress, and the deformation becomes a plane strain. 
The elastic response of the enclosed area is given by 
Eei = (Ti — va2 and Ee2 = a2 — i'CTi, where ei, 62 and cti, 
(T2 are the plane strain and normal stress components, 
respectively. It immediately follows that the relative 
change in the square area equals ei-|-e2 = —{\ — v)E~^2p 
where p = —(0-1 + cr2)/2 is the two-dimensional pressure. 
The enclosed area deforms under this pressure into an- 
other square so that a constant gap 5g is created along 
the loop. Using that p = f{x)dx, the size of this virtual 
gap becomes 



(5g = (1 - v)E-^xf{x)dx 



(10) 



The sign of 5g is so that a positive value represents a void 
space, whereas a negative 5g corresponds to an expansion 
of the inner square. 

The contribution from this elementary deformation to 
the strain of the whole body is obtained by recombining 
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FIG. 3. (a) Partial effect of the volume force in a square 
cylinder. By the shown transformation the elastic problem 
can be treated analytically using superposition, see text, (b) 
Schematic picture of the deformation (strongly amplified) of 
a square cylinder in the case of (left) incomplete and (right) 
complete flux penetration. 

("gluing together") the parts on each side of the virtual 
gap. The resulting strained state is equivalent to hav- 
ing a pair of edge dislocations in each corner of the loop. 
As each pair can be replaced by one dislocation directed 
along the diagonal, we conclude that the original gap 
space can be replaced by two slits of width 8G — 5g^/2 
forming a central cross. By this transformation the in- 
tegral effect of the entire force field can be determined 
analytically as follows. Let y denote the distance from 
the center to a point on a diagonal. At this point the 
total gap, G, is an accumulation of the gaps 5G created 
from y and out to the corner y = w\l. The total result 
is therefore given by 



G(y)^^l-. 



w 



V2 



E^^0 Jy/V2 



B{x)B'{x)xdx 



(11) 



The full deformation is now obtained by "gluing" the vir- 
tually disrupted material together. Since G'(jj) ex Bjcy 
it follows that G{y) is not a linear function, thus prov- 
ing that the pinning-induced magnetostriction does not 
conserve shape in the square case. 
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FIG. 4. Gap shape at various stages of flux penetration, 
(a) Increasing field {Ba indicated) and (b) decreasing field for 
Ba/Bp = 0.8, 0.5, 0.2, (thick line), -0.2, -0.5, -0.8 -1.1. 



G(i/) with an oscillatory shape, thus giving a mixed type 
of deformation. As the reverse field increases in magni- 
tude the shape again becomes concave. For Ba < —Bp 
the sample shape follows the same development as during 
the field increase. 

In conclusion, we have investigated and compared the 
pinning-induced magnetostriction in two important ge- 
ometries. In contrast to the shape-conserving circular 
case, which could be solved exactly, the existence of cur- 
rent discontinuity(d)-lines in the square case complicates 
the behavior by producing shape distortion. We have 
predicted the appearance of convex, concave and even 
more complicated types of deformations where the dis- 
tortion is of the same order of magnitude as the over- 
all striction. Since in most cases HTSC crystals have a 
rectangular shape, which always lead to d-lines in the 
critical-state, the distortion could easily lead to misin- 
terpretation of dilatometric measurements. The problem 
should be minimized by (i) having a high aspect ratio 
rectangular shape so that the central part of the sample 
is described approximately by a plane strain deforma- 
tion, and (ii) designing the displacive sensor with a local 
contact to the long side of the rectangle in order to avoid 
distorted corner regions. Full analysis of the proper cor- 
rections will clearly require numerical work. 

The authors acknowledge discussions with Yu. 
Galperin, and financial support from The Research 
Council of Norway. 



The eq.(|l|) allows the gap to be calculated. In the 
Bean model along the increasing field branch one finds 
for G = G EiJ.o{l — u)~^B~'^ the expression 



G{v) 
wV2 
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I) + 2-2(^)^ (12) 
w^/2 



in the penetrated region. In the central Meissner area the 
gap is constant. The behavior is shown in Fig. 4 a, and 
representative examples of the resulting concave defor- 
mations are illustrated qualitatively in Fig^ b. For the 
fully penetrated state one sees from eq.(|2|) that G{y) 
grows linearly with B^, while at the same time G — Q 
at the corners. This produces a deformation where the 
concave character grows steadily with the applied field. 

The behavior of G{y) for descending applied field is 
shown in Fig. 4 b. The deformation along this branch 
is evidently more diverse even when we omit the stage 
following immediately after field reversal. The shown 
graphs represent only stages where jc is reversed through- 
out the sample. The profiles with negative G gives rise to 
convex deformations. The convexity persists all the way 
down to the remanent state, where the distortion near 
the corners is vanished. Then, as the field is reversed 
there is an interval where the behavior is governed by a 
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